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ABSTRACT 

Context. The extremely relativistic Fe lines, detected from some active galactic nuclei (AGN), indicate that generation and repro- 
cessing of the X-ray emission takes place in the immediate vicinity of the event horizon. Recently, general relativistic (GR) effects, 
in particular light bending which is very strong in that region, have been considered to be the cause of complex variability patterns 
observed in these AGNs. 

Aims. We study in detail the GR models of the X-ray spectral variability for various geometries of the X-ray source and with various 
relativistic effects being the dominant cause of spectral variability. The predicted properties are compared with the observational data 
of the Seyfert 1 galaxy MCG-6-30-15, which is currently the best studied AGN with signatures of strong gravity effects. 
Methods. We focus on modeling the root mean square (RMS) spectra. We compute the RMS spectra for the GR models using a Monte 
Carlo method, and compare them with the RMS spectra from the Suzaku observations of MCG-6-30-15 on January 2006. 
Results. The data disfavor models with the X-ray source (1) moving vertically on the symmetry axis or (2) corotating with the disc 
and changing height not far above the disc surface. The most likely explanation for the observed fractional variability is given by the 
model involving the X-ray source located at a very small, varying distance from a rapidly rotating black hole. This model predicts 
some enhanced variations in the red wing of the Fe line, which are not seen in the Suzaku observations. However, the enhanced 
variability of the red wing, while ruled out by the Suzaku data, is consistent with an excess RMS variability, between 5 and 6 keV, 
reported for some previous ASCA and XMM observations. We speculate that the presence or lack of such a feature is related to the 
change of the ionization state of the innermost part of the disc, however, investigation of such effects is currently not possible in our 
model (where a neutral disc is assumed). If the model, completed by a description of ionization effects, proves to be fully consistent 
with the observational data, it will provide a strong indication that the central black hole in MCG-6-30-15 rotates rapidly, supporting 
similar conclusions derived from the Fe line profile. 

Key words, galaxies: individual: MCG-6-30-15 - galaxies: active - galaxies: Seyfert - X-rays: galaxies 

1. Introduction The time-resolved X-ray spectra in MCG-6-30-15 indicate 

that the best explanation of its spectral variability pattern is given 

nr^-KT f , f ^ , ■ ui V • • by a phenomenoloeical, two-component model consisting of (i) 

AGNs are powerful sources of strongly variable X-ray emission, , . , , ■ , • , , ^. ' 

.,•11 ■ 1 » » 1 ui 1 u 1 T-u • a highly variable power-law continuum (referred to as primary 

most likely originating close to a central black hole. The on- ■ ■ \ a f \ u i ■ v.^ a ^ a 



■ gin of the X-ray variability is poorly understood, however, some 



emission) and (ii) a much less variable, and uncorrelated, com- 



S3 - J,., , A • J » u j' *i ^\ A ^ ponent with a spectrum characteristic of radiation reflected from 

of the observed variations are supposed to be directly related to ^ . ,. , , ■ n \ — rrr ; „ „ . . — i 

■ ■ V . J »u • • • • • u» the inner disc (e.g. Fabian & VaugharJ 2003UVaughan & FabiarJ 

strong gravity effects and their investigation may give insight ■^r--— , r; ^ " , — — = — . *\ ^. f 

• . .. f .. » • • »u • • f 7u 2004, Larsson et al. 2007; an alternative explanation, proposed 

into properties of the space-time metric in the vicinity of the , , ' „ , . , . „ ,. ■ . r, ™ , 



event horizon by Miller et al. 2008, is briefly discussed in Sect. \5M- The lat- 
ter component includes both the Fe Ka line and the associated 

The bright Seyfert 1 galaxy MCG-6-30-15 (z = 0.00775) is Compton reflection hump, peaking ai'ound 30 keV, produced by 

the key object for such studies. Its X-ray spectrum shows an ex- down-scattering of higher energy photons. The varying relative 

ti-emely distorted Fe Ka line, indicating that X-ray repr ocessing contributions of these spectral components can explain the ob- 

takes place very clos e to the central black hole (e.g. Fabian et all served hardening of the X-ray spectrum at lower fluxes as well 

l255llMiniuttiirn]l2057h . Furthermore, crucially for investiga- reduced variations of flie blue peak of die Fe Ka hne. 

tion of strong gravity in MCG-6-30-15, radiation reflected from The apparent lack of connection between the two compo- 

the inner disc is not contaminated by radiation reflected from dis- nents indicates that some complex physical mechanism should 

tant matter (see Sect. 14. 2t . The disc-line model for the detected be involved, as the components are supposed to originate 



line profile requires that roughly half of the observed Fe photons together in the central region and a strict correlation be- 
come from the innermost parts of the accretion disc, within 6 tween them would be ex pected in a simple reflection picture, 
gravitational radii (these photons form the red wing of the line, iFabian & Vaughaiil (1200 3') first argued that this disconnectedness 
below 6 keV) and the remaining from a slightly more distant re- may be explained by relativistic effects, in particular by light 
gion, at a few tens of gravitational radii (forming the blue peak, bending and focusing of the primary emission towards the ac- 
observed between 6 and 7 keV). cretion disc. In this scenario, elaborated under some specific as- 
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sumptions by iMiniutti & FabianI ( |2004 MF04), the change of 
the height of the primary X-ray source is the cause of spectral 
variations. 

However, a systematic inves tigation of this class of models 
bv lNiedzwiecki & ZvckH ( l2008l NZ08) questioned their ability 
to explain the observed effects. In particular, MF04 assess the 
reduced variability of the reflected component based on the be- 
havior of the total flu x of fluorescent photons received by a dis- 
tant observer. iNZOSi point out that changes of the height of the 
source induce substantial changes in the line profile, with sig- 
nificant variations of the fluxes in the blue peak and in the red 
wing. These balance each other for some range of parameters, 
reducing variations of the total flux, however, this property is 
not sufficient to explain the observed data. Then, the model fails 
to reproduce the key property which motivated its development, 
i.e. reduced variation of the blue peak of the Fe line. 

Interestingly, NZ08 find an alternate model, which can ex- 
plain reduced variations of the blue peak. The model propose d 
in NZ08 follows the generic idea of iFabian & VaughanI (l2003h . 
involving a scenario where the change of position of the primary 
source leads to large changes of directly observed emission and 
smaller changes of the amount received by the disc. However, 
the model assumes changes of the radial distance rather than the 
height and its properties rely on a qualitatively distinct effect, 
namely bending to the equatorial plane of the Kerr space-time. 
The model proposed in NZ08 predicts, however, some enhanced 
variations of the red wing of the Fe line, which may challenge 
its applicability. 

At present, the time-resolved spectra for sufficiently short 
(~ 10 ksec) time bins have poor quality, which does not allow us 
to study such fine details of spectral variability. Thus, less direct 
tools are used. In particular, the RMS spectra - sh owing the frac- 
tional variabi lity as a f unction of energy (see e.g.'E delson et al.l 
[2002; Marko witz eta0 2003: Vaughan et al. 2003b) - aUow us to 
investigate spectral variability in a model-independent manner 

In MCG-6-30-15, the general trend is that the RMS spec- 
trum decreases with the increase of energy, moreover, it has a 
pronounced, broad depression around 6.5 keV. Such a shape of 
the RMS spectrum has been confirmed by the observational data 
from ASCA (Matsumoto et al. 2003), XMM (Vaughan & Fabian 
2004) , j^Xrg (Markowitz et al. 2003) and Suzaku (.Miniutti et ah 



2007h . IVaugh an & Fabian (2004 show that fliis form of the 



RMS spectrum may be reproduced in a model with a variable 
power-law and a constant reflection. Again, results of NZ()8| 
challenge the physical motivation for such a phenomenologi- 
cal description in terms of GR models, as these typically pre- 
dict some changes of the spectral shape of reflected component. 
Changes in the Fe line profile have been studied, including their 
impact on the RMS variability, in NZ08. However, that paper ne- 
glected the Compton reflected component and its results cannot 

be directly compared with the observed data. 

In this paper, we extend the model developed in lNZOSl bv in- 
cluding the Compton reflected component, and we analyze its 
predictions for the fractional variability amplitude. Then, we 
compare the simulated RMS spectra with the spectra derived 
from three Suzaku observations of MCG-6-30-15 performed in 
January 2006. 



2. GR models of the X-ray spectral variability 

2.1. The GR models 

We consider an accretion disc, surrounding a Kerr black hole, ir- 
radiated by X-rays emitted from an isotropic point source (here- 



after referred to as the source of primary emission). The black 
hole is characterised by its mass, M, and angular momentum, J. 
We use the Boyer-Lindquist coordinate system x' - (t, R, 6, (p) 
and the following dimensionless parameters 



R 

rI 



J 



cRgM 



df 



(1) 



where Rg = GM/c is the gravitational radius and f = ct/Rg. 

We assume that a geometrically thin, neutral, optically thick 
disc is located in the equatorial plane of the Kerr geometry. For 
distances greater than the radius of the marginally stable circular 
orbit, r„is, we assume circular m otion of matter formi ng the disc, 
with Keplerian angular velocity dBardeen et al.ll 19721) . 
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(2) 



We take into account reflection from matter which free falls 
within Tms, however, this effect is negligible for the large val- 
ues of a considered in this paper. We assume the outer radius 
of the disc rout - 1000. Negligibly small dips, related to this fi- 
nite outer extent, occur in the simulated RMS spectra around 6.4 
keV. 

Inclination of the fine of sight to the rotation axis of the black 
hole is given by flobs- Our discussion of properties predicted by 
GR models, in Sect. [3] focuses on values of the inclination an- 
gle between Oobs -25° and 45°, which is relevant to modeling 
MCG-6-30-15. The location of the primary source is described 
by its Boyer-Lindquist coordinates, rj and 6^, or by = rj cos 0^ 
and ps = '"s sin 0^ . Note that the azimuthal position of the source 
is not relevant to our results, as we consider only spectra aver- 
aged over a complete orbit of the source, see Sect. 12.41 

We neglect transversal or radial motion of the source. We 
take into account the azimuthal motion with the generic assump- 
tion that the X-ray source (if displaced from the symmetry axis) 
corotates with the disc. We note, however, that the description of 
such a corotating source is ambiguous at larger and the choice 
of a specific assumption may significantly affect properties of 
the model. Plausible assumptions include relating the angular 
velocity of the source, Qs('"s, ^s), to the Keplerian angular veloc- 
ity in the disc plane at or at ps. In general, the former, with 
QK(rs), yields an azimuthal velocity of the source, vf, slightly 
smaller than the Keplerian velocity of the disc, (^ 0.5c in the 
innermost disc), while the latter, with Qk(Ps), yields vf signifi- 
cantly larger than v^; note that the azimuthal velocities, vf and 

v^, are defined wi t h res pect to the locally non-rotating frame, 
see i Bardeen et al.l ( 119721) . In particular, for hs/r^ - 1, using 
Q.S - flK('"s) we get vf ^ 0.4c, while Q.^ - Ok(Ps) yields vf 
exceeding 0.6c and 0.8c at rs = 3 and - 2, respectively. The 
above issue is particularly relevant in model C, defined below, 
where Qj = Q.K{Pf,) is assumed (following MF04). In turn, our 
model S (see below) assumes Qj - Qk('"s), although at small 
/js/^s, considered in this model, the difference between the two 
approaches is negligible. 

The relevant relativistic effects are clearly represented in the 
following three models. 

Model S involves a source located close to the disc sur- 
face and rotating with the Keplerian velocity, Q.K{fs)- We con- 
sider very small heights of the source above the disc surface, 
/js Ts, for which rigid coupling (and corotation) of the source 
is the most likely configuration; moreover, the description of az- 
imuthal motion is not subject to the ambiguities noted above. 
Variability effects in model S result from varying radial distance. 
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Ts, at a constant polar angle 6^ - Note that the change of the radial 
position leads to slight changes of the height, hs, so the source 
motion is conical rather than plane parallel. 

Most results presented for model S in this paper are derived 
for a fixed polar angle, 6^ = 1.5 rad, yielding hs/fs - 0.07; in 
Sect. [33] we show results for other values of 6^. 

Two different physical effects can be studied in the following 
two regimes of model S . At large distances, Ts > 6, the bulk of 
the reflected radiation arises from a small spot under the source; 
moreover, the gravitational redshift is weak and the variability 
properties are dominated by local Doppler distortions, see Sect. 
13.21 For < 4 and high values of a (> 0.9), a qualitatively dis- 
tinct property of the Kerr metric (namely the bending of photon 
trajectories to the equatorial plane) results in both strong and 
approximately constant illumination of the regi ons of the disc 
where the blue peak of the Fe line is formed, see lNZOSl The lat- 
ter case, i.e. model S with large a and small r^, will be referred 
to as model S^'^. 

Model A involves a static primary source located on the 
symmetry axis. Such a scenario, often referred to as a lamp- 
post model, has been considered in a number of papers, e.g. 
[Martocchia & Matt ( 1996), Petrucci & Henri ( 1997). Variability 
in model A is induced by the changes of the source height, h^, 
above the disc. This model allows us to investigate properties re- 
lated to bending of photon trajectories toward the center, a dom- 
inant effect for such a location of the source. 

Our model C corresponds to a variability model with a 
cylindr ical-like motion of the primary source, developed by 
|MF04| (note that their original model assumed a small value of 
Tout = 100, w hich influenced some of their results, see discus- 
sion in|NZ08). The model (often referred to as the light-bending 
model) involves a source changing its height above the disc, sim- 
ilarly to the lamp-post model, but displaced from the axis, with 
the constant projected radial distance p,, and rotating around it. 
Furthermore, this model assumes that at each the source has 
the same angular velocity, Qk(Ps)- As noted above, the last as- 
sumption yields vf significantly exceeding the Keplerian veloc- 
ity of the disc. Moreover, the dependence of vf on /j, is non- 
monotonic, which is crucial for the variability predicted by this 
model, see iNZOSl Then, model C allows us to study the impact 
of kinematic effects, which are complementary to the gravity ef- 
fects underlying properties of models A and S . All results for 
model C shown in this paper are derived with ps = 2. 

In models A and C we consider the maximum rotation of a 
black hole, with a - 0.998. For model - which appears to 
be the most feasible to explain the observed data - we investigate 
its dependence on the value of a. 

2.2. Intrinsic luminosity 

In models A and C we assume the same intrinsic luminosity of 
the primary source at all h^. Model S , in Sect. 13.21 also assumes 
an identical luminosity at various r,. For model S^^ we focus, 
however, on a (physically motivated) scenario where the radial 
profile of intrinsic luminosity, L^(r s), follows the diss ipation 
rate per unit area in a Keplerian disc iPage & Thornel[T97 4). All 
computations for models A and C, and for model S or 5 '^^ in 
Sects. 13.2143.41 are performed assuming a constant (in time) lu- 
minosity at any given location. For model 5^^, we consider two 
further modifications: 

(i) A power-law modulation of the radial luminosity profile, rf x 



(ii) As a step toward more realistic modeling, we take into ac- 
count variations (in time) of intrinsic luminosity. In most cases 
we assume a Gaussian distribution of luminosities, with the de- 
viation ctl, but we note that properties of the model depend on 
the assumed distribution function; details are given in Sect. 13.51 

2.3. Monte Carlo computation of the RMS spectra 

For each model, we compute the observed spectra for a number 
of primary source positions, using the following grids of or /is. 
In models A and C the heights are linearly spaced with Ah^ = 1 
at /is > 3 (so /is - 3, 4, 5, ...) and A/is = 0.2 at /is < 3. In model 
S the radial distances are linearly spaced with Ars = 1 ; however, 
in the regime of model S^^ (i.e. at rs < 4) we use Ars =0.1. 

For a given position of the primary source, we use a Monte 
Carlo method, involving a fully general relativistic treatment of 
photon transfer in the Kerr space-time, to find spectra, observed 
by a distant observer, of both the primary emission and the re- 
flected radiation. A large number of photons are generated from 
the primary source (for each position we typically trace 10^ pho- 
tons) with isotropic distribution of initial directions in the source 
rest frame and a power-law distribution of photon energies, char- 
acterised by a photon spectral index, F. Results presented in this 
paper correspond to F = 2. We have not found any noticeable 
difference in the RMS spectra computed for various spectral in- 
dices in the range relevant to MCG-6-30-15, F = 2 - 2.3 (in- 
dicate d by modeling of time-resolved spectra, e.g. lLarsson et al.l 
l2007l) . 

For each photon, equations of motion are solved to find 
whether the photon crosses the event horizon, hits the disc 
surface or escapes directly to the distant observer For pho- 
tons hitting the disc, we perform a Monte Carlo simulation of 
Compton reflection. Photons transferring in the disc are sub- 
ject to consecutive Compton scatte ring events competing wit h 
absorption. We use abundances of lAnders & Grevessd ( Il989l) . 
For an isotropic illumination of the disc surface, our code repro- 
duces, in the disc rest f rame, the reflection spectru m described 
by the pexrav model (Magd ziarz & Zdziarski|[T995h . Our simu- 
lation of Compton reflection allows us to take into account an 
incidence-angle dependent irradiation of the disc. 

For a photon hitting the disc surface with energy > 7.1 keV, 
we generate an iron Ka photon, with energy 6.4 keV, emerg- 
ing from the disc. The relative weight of the Fe photon is re- 
lated to the initial energy and direction of th e incident photon 
by the quasi-analytic formu la [eqs. (4-6)] from lGeorge & FabianI 
(1991"). Similarly to 'NZ08, we modify the original formula by 
multiplying it by a factor of 1.3 to account for elemental abun- 
dances, consistent with those assumed for Compton reflection. 
We assume a limb-darkened emission of Fe Ka photons, with 
intensity /(//em) 1 + 2.06pem, where yUgm = cos 0em and ^em 
is the polar emission angle in the disc rest frame. As shown in 
INZ081 for a limb-brightened emission, only small, systematic 
c hanges - related to the strength of the blue peak (see section 3.7 
in lNZ08h - occur in the RMS spectra. We solve equations of mo- 
tion in the Kerr metric for both the Fe Ka and Compton reflected 
photons. We take into account reflection of photons that return 
to the disc, following the previous reflection. 

We compute the observed spectra for a number, A^, of various 
positions of the primary source and then we determine the RMS 
spectrum according to the following definition 



RMS(£) = 



1 



< m > \ 



[/;■(£)- < f{E) >]2 

N-\ 



(3) 
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where fi{E) is the photon flux in the energy band, E, correspond- 
ing to /-th position of the source and < f(E) > is the average 
photon flux in this band for all positions. 

In our basic analysis, in Sect. 13.1143.31 we assume that each 
spectrum, //, corresponds to a specific position of a point source 
and the RMS spectra are derived taking into account all posi- 
tions within a range specified by its lower and upper ends, rmin 
and r,nax or /i^in and /imax, respectively; this typically involves 
summation over ~ 10 positions of the primary source. 

In such a case (i.e. for spectra, fi, from a single position), the 
impact of relativistic effects on spectral variability is maximized 
(see Sect. 13. 4t . However, such a large magnitude of relativistic 
effects can be observationally revealed only under very specific 
conditions. Namely, the light curve should be sampled in time- 
bins with a size comparable to (or shorter than) the time-scale 
of the change of position of a primary source. Furthermore, ap- 
proximation of the hard X-ray emitting region by a single, point- 
like source may be inadequate. Particularly in model S , multiple 
flares may occur simultaneously at various sites on the disc. 

To estimate predictions of GR models for more complex ar- 
rangements of the primary source, or larger time-bin sizes, we 
study the RMS spectra for the model with the energy spectra, 
fi, formed by contributions from several random positions of the 
source. We make a simplifying assumption that each basic spec- 
trum, fi, contains a contribution from the same number, nx, of 
locations. Each location is randomly generated with a continu- 
ous probability distribution; then, or from the grid of a given 
model, nearest to the generated value, is taken into account. In 
models A and C, /is is generated with a uniform probability dis- 
tribution between /imin and /imax- In model S^'^ we use the prob- 
ability density in the form of a power-law, P(rs) fi- The RMS 
spectra in these more complex cases are computed with = 30. 

Summarizing the above, model S which is studied most 
thoroughly here, follows the kinematic assumptions of model S 
and, unless otherwise specified, assumes a constant intrinsic lu- 
minosity given by L^(rs). Modifications of these basic assump- 
tions are parametrized by: /3 - for the change of the radial lumi- 
nosity profile; <tl - for variations of intrinsic luminosity at (each 
given) Ts; and 6 - for the radial probability distribution. 

2.4. Azimuthal averaging; time deiays 

We ignore time delays between the primary and reprocessed ra- 
diation, related to light-travel times between the source and var- 
ious parts of the disc; moreover, we determine the RMS spectra 
using the spectra, fi, averaged over azimuthal angle. We empha- 
size that both of these issues are not important to our analysis, 
which focuses on effects observed in MCG-6-30-15 (harbouring 
a black hole with a mass, M ^ 3 x 10''Mo; McHardy et al. 2005.) 
on a time-scale of 10-100 ksec. 

Regarding the time delays, we note that - in most cases - 
almost all the reprocessed radiation originates from r < 30 
and the delays are smaller than 1 ksec. Only in model C with 
/is > 10, irradiation of the disc at r ^ 100 cannot be ne- 
glected, for which time delays of a few ksec occur. These would 
be relevant for modeling some effects studied for short time- 
scales, e.g. in poin t-to-point RMS spectra with 1 ksec bins (cf. 
IVaughan & Fabianll2004b . The time bins used to derive the RMS 
spectra in this paper are an order of magnitude larger and any 
variations on the time-scale of a few ksec would be grossly un- 
dersampled. 

Similarly, in most cases, the orbital period of the source 
does not exceed 1 ksec, which justifies the azimuthally-averaged 
treatment. The only exception involves model S with large 
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Fig. 1. Panel (a) shows the RMS spectra for model A with /zj 
changing between 1.8 and 10. Panel (c) shows the RMS spectra 
for model C with /is changing between 4 and 20. In both pan- 
els, the spectra from top to bottom are for flobs = 25° (solid, 
green online), 30° (dashed, black; shifted vertically by -0.02, 
for clarity), 40° (dotted, red; shifted vertically by -0.04) and 45° 
(dot-dashed, blue; shifted vertically by -0.06). Panels (b) and (d) 
show the RMS spectra for various spectral components in model 
A and C, respectively, with flobs = 30° and other parameters indi- 
cated for panels (a) and (c). The curves, from bottom to top, are 
for Compton reflected hump (dotted, blue; note that it does not 
include the Fe line), primary 4- Compton reflection (solid, red) 
and primary (power-law) component (dot-dashed, magenta). The 
dashed (black) curve is for the total emission, i.e. primary con- 
tinuum + Compton reflection + Fe line. 



(^ 10). We neglect the detailed study of this regime, as it ap- 
pears not relevant in modeling the observed data. 



3. Results 

Figures [1] [3]|4] and |5] show the RMS spectra for models A, C 
and S , derived under the assumption that each of the individual 
spectra, fi(E), contributing to the RMS spectrum corresponds to 
a fixed, single position of a primary (point) source. The RMS 
spectra are computed for a constant intrinsic luminosity of the 
primary source. Figures|6]and|2]show the RMS spectra for more 
complex scenarios, discussed in Sects. [3741 and l3. 51 



3.1. Ctianges of hs 

As shown in Figs.[Tla)(c) and[3] models assuming a vertical mo- 
tion of the primary source do not predict a depression around 6.5 
keV, similar to that revealed in MCG-6-30-15 (cf. Fig.[8]below). 
Moreover, an opposite property, i.e. increase of fractional vari- 
ability at this energy, is typically predicted by these models in 
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Fig. 2. Changes of the Fe line profile in model A (panels a and c) 
and model C (panels b and d). The top panels are for 0obs = 25°; 
the bottom panels are for ^obs - 40°. In all panels, the profiles 
from bottom to top are for /is = 4, 6, 8 and 10. Note that, for 
^obs = 40°, rather subtle differences around the maximum of the 
blue peaks lead to opposite properties, i.e. a sharp excess or a 
sharp drop, around 6.9 keV in the RMS spectra for model A and 
C, respectively, see Fig.[T] The units are arbitrary but the same in 
all panels; profiles in panels (a) and (b) are rescaled by a factor 
of 0.55. 
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Fig. 3. RMS spectra in two regimes where the fractional variabil- 
ity does not decline with increasing energy. The bottom (solid, 
black online) is for model C with changing between 0.2 and 
1.2. Only in this range of parameters is the reflected radiation 
more variable than the primary emission and a pronounced fea- 
ture, related to the varying shape of the Fe edge, occurs in the 
RMS spectrum. The upper (dashed, red) is for model S with rj 
changing between 6 and 12. Both spectra are for 0obs = 30°. 



most of the parameter space. Before discussing the specific mod- 
els, we comment on the generic impact of various spectral com- 
ponents on formation of the fractional variability spectrum, illus- 
trated in Fig. [nb)(d). The narrow excesses in the RMS spectra 
result from changes of the Fe line profile. These are investigated 
in detail in NZ08; here we briefly discuss some basic properties 
crucial for this study. We point out that relativistic distortions of 
the Compton hump are not as crucial for the fractional variability 
as distortions of discrete spectral features. This property results 
from a continuous energy distribution of the hump and, in most 



cases, leads to decrease of the RMS spectrum, consistent with 
observations. 

The primary continuum is the most variable component, re- 
garding variations of the directly observed flux, but it does not 
change the spectral shape. Then, its variability is represented by 
a flat RMS spectrum with a large amplitude. In turn, changes of 
the amount of primary radiation illuminating the disc - resulting 
from the change of /zj - are smaller (as a result of light bending). 
However, the distance of the most intensely illuminated area of 
the disc changes with /is. As a result, discrete spectral features 
in the reflected radiation, in particular the Fe Ka line and edge, 
are subject to varying relativistic distortions, yielding an excess 
fractional variability around their rest energies. An excess, re- 
lated to the variability of the blue peak of the Fe Ka line, is seen 
in most of the spectra in Fig. [na)(c) (except for model C with 
large ^obs, see below). As we noted before, the excess variability 
caused by the Fe Ka line occurs even when the total flux in the 
line does not change significantly. 

Similarly, changes of the shape of the Fe Ka edge result 
in a pronounced excess, around 7 keV, in the RMS spectra for 
the Compton reflected component, see the bottom curves in Fig. 
[nb)(d). However, variability at this energy is dominated by the 
primary continuum and the Fe line, and the excess (due to the Fe 
edge) does not yield any significant signature in the RMS spec- 
trum for the total observed radiation (except for model C with 
small hfi, see below). Crucially, though, for interpretation of the 
observed data, the distortions of the Fe edge yield a smooth RMS 
spectrum, see the middle solid curves, for the sum of the primary 
continuum and the Compton reflected radiation, in Fig.[nb)(d), 
unlike the case of reflection from a static slab, for which a sharp 
drop occurs at 7.1 keV (cf Fig.[8]below). 

The Compton hump, although subject to the same relativistic 
distortions as the Fe line, is less significantly affected by their 
changing amount, owing to the continuous spectral distribution 
of the hump. As a result, the amplitude of the RMS spectrum for 
the Compton reflected component is typically much smaller than 
for the primary emission, compare the top and bottom curves 
in Figs. [nb)(d). As the contribution of the Compton reflected 
radiation to the total spectrum increases between 2 and 30 keV, 
the RMS spectrum for the total emission decreases, reaching a 
minimum around 30 keV where the reflection hump peaks. 

3.1 .1 . Model A: light-bending (to the center) 

Trajectories of photons emitted from a source located close to 
the symmetry axis are subject to a simple bending to the center, 
independent of the value of a (contrary to bending to the equa- 
torial plane, discussed in Sect. |3.3l l. Models involving a verti- 
cally moving source close to the axis predict strong variations of 
the observed flux. They result from changing amounts of purely 
GR effects, i.e. the light bending and gravitational time delay. 
However, these models robustly predict an enhanced variability 
around the Fe Ka energy, see Fig. \lla). This property, incon- 
sistent with observations, property results from the fact that a 
source located at a larger hs illuminates a more extended area 
of the disc and, therefore, the increase of /is gives rise to rather 
significant strengthening of the blue peak (see Fig.|2^). The vary- 
ing strength of the peak yields a pronounced excess in the RMS 
spectrum between 6 and 7 keV. 

Dependence on 0obs, seen in Fig.[Tfa), results from gravita- 
tional focusing and Doppler coUimation, toward observers with 
high inclination angles, of radiation reflected from the inner- 
most parts of the disc (which receive most of irradiating flux 
from a source located at small hs). A stronger contribution of 
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the Compton reflected radiation for larger 0obs gives rise to a 
stronger decrease of the RMS spectra with increasing energy. 
The focusing toward higher 0obs results also in weaker variations 
of the blue peak (compare Figs.|2ja) and|2lb)) and, thus, in a less 
pronounced RMS excess around 6.4 keV. 



3.1.2. Model C: azimuthal motion 

The major difference between models A and C results from rota- 
tion of the primary source around the symmetry axis in the latter 
model. Then, in model C, changes in the fluxes of the directly 
observed and the irradiating radiation, induced by the change of 
/is, result from the combination of light bending and Doppler 
beaming related to this motion (we recall that in both models 
A and C the intrinsic luminosity remains constant). The specific 
kinematic assumptions of model C yield a non-monotonic de- 
pendence of azimuthal velocity on h^, leading to different proper- 
ties of the model in two ranges of (see NZ08 for details). The 
critical height, approximately limiting these two ranges, is the 
one at which the azimuthal velocity achieves a maximum value 
(for Ps - 2 this maximum occurs at /zj = 2). At small /is (< 2 for 
Ps - 2), variations of the directly observed primary continuum 
are strongly reduced and the RMS spectrum is shaped by spectral 
changes of the Compton hump; note that the reduced variations 
of the primary flux are observed only at smaller 0obs (i-C- these 
considered here) and result from more efficient Doppler beam- 
ing to larger ^obs, due to the increase of vf , which approximately 
balances less efficient light bending with the increase of h^. The 
predicted increase of the RMS variability for increasing ener- 
gies, accompanied by excess variability around 6 keV (caused 
by the changing Fe Ka edge), see Fig. [3] is clearly not consis- 
tent with observed data. 

At larger heights, the trends in both components are reversed. 
The primary continuum is now the most variable component, 
yielding the decrease of the RMS spectrum towards higher en- 
ergies. Effects related to changes of the line profile strongly de- 
pend on ^obs- The crucial property resulting from the Doppler 
beaming of primary emission is that the reflected radiation orig- 
inates from more extended regions of the disc (with r > h^, as 
opposed to model A) and hence the Fe line is somewhat less vari- 
able than in model A (see Fig. |2]i. This is reflected in the RMS 
spectrum by the sharp drop, with the width of a few hundred e V, 
occurring for 0ob,s > 30° at the energy of the maximum of the 
blue peak. Regardless of the value of flobs, model C predicts a 
pronounced depression around 5.5 keV. 



3.2. Model S, large r^: local Doppler effects 

The upper curve in Fig.[3]shows the RMS spectrum for model S 
with large Ts (> 6). Radiation reprocessed from the emission of 
a source located close to the disc surface at large rs originates in 
bulk from a small spot below the source. In such a case, both the 
directly observed primary emission and the reflected radiation 
are subject to the same relativistic distortions and, therefore, the 
relative normalization of the latter component is the same as for 
a static slab illuminated in flat space-time. The departure of the 
RMS spectrum from a flat form results from varying Doppler 
distortion of the Fe line, corresponding to changes of rs. These 
changes yield a strong excess around 6.4 keV due to the change 
of the energy of the maximum of the blue peak. An additional, 
smaller excess related to the varying extent of the Une occurs at 
lower energies. 
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Fig. 4. Panel (a) shows the RMS spectra for model S ^'^ with dif- 
ferent values of spin parameter a = 0.9 (dotted, blue online), 
0.95 (dashed, magenta), 0.98 (solid, red) and 0.998 (dot-dashed, 
black), from top to bottom. For all curves, rs varies between 
'"mill - '"ms(fl) and rmax - 3.2, and 0obs = 30°. Panel (b) shows the 
RMS spectra for various spectral components, as in Fig.[nb)(d), 
in the model with a = 0.98. 



Note, however, that such a strong excess corresponds only to 
a simplified scenario with nx - I (see Sect. |3.4| l. 

3.3. Model S^^ (small in the Kerr metric): bending to the 
equatorial plane 

Figures|4]and|5]show the RMS spectra for model S^'^ with small, 
varying rs and large a. Such parameters define a unique regime 
with the reduced variability of the blue peak of the line occurring 
independently of other specific assumptions. In this scenario, the 
reflected radiation has two components with different variability 
behavior. The first one arises locally from a strongly irradiated 
spot (or - after azimuthal averaging - a narrow ring) under the 
primary source. In principle, this component follows the behav- 
ior described in Sect. l3.2l but - in addition - it is subject to strong 
(and varying with rs) gravitational redshift. The variable redshift 
of the Fe line gives rise to a pronounced excess in the RMS spec- 
trum between 4 and 6 keV. The second component arises from 
a slightly more distant region, at ~ lORg, which is strongly il- 
luminated due to bending to the equatorial plane, an effect sig- 
nificantly affecting radiation emitted from rj < 4 for a > 0.9. 
It is the second component which yields the unique properties 
of this model, through the combination of the following two ef- 
fects: (1) the observed blue peak is formed mostly by photons 
emitted at ~ lOR^; and (2) the flux illuminating that site remains 
approximately constant while rs changes (see NZ08 for details). 
The resulting reduction of variations of the blue peak yields a 
pronounced depression in the RMS spectrum between 6 and 7 
keV. As the reduction effect fully relies on the properties of the 
Kerr metric, it is obviously stronger for larger a, see Fig. Ufa). 

Similarly to model A and model C with hs > 2, the Compton 
component is less variable than the primary continuum, see Fig. 
|4f b), and as a result the RMS spectrum decreases with increasing 
energy. 

The most significant change in the RMS spectrum, corre- 
sponding to the change of flobs, is related to the shape and loca- 
tion of the excess between 4 and 6 keV (produced by variations 
of the red wing), see Fig.|3a). Similarly, the change of 9^ is most 
significantly reflected in the shape of the excess. For smaller 9^, 
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Fig. 5. Dependence of the RMS spectra in model S^^ on: 0obs 
(a); 6*5, or equivalently hjr^ (b); r^i^ and r^ax (c); /? (d). All 
spectra are for a = 0.998; panels (b-d) are for 0obs = 30°; pan- 
els (a),(b) and (d) are for r^in = 2 and r^ax = 3. Panel (a): 
^obs = 25° (dot-dashed, magenta online; shifted vertically by 
0.06), 30° (dotted, blue; shifted by 0.04), 40° (solid, red; shifted 
by 0.02) and 45° (dashed, black). Panel (b): the spectra from top 
to bottom are for h^/r^ = 0.07 (solid, red; shifted vertically by 
0.04), hjr, = 0.14 (dashed, black; shifted vertically by 0.02) 
and /zs/^s = 0.21 (dotted, blue). Panel (c): the spectra from top 
to bottom are for rs varying between r^in = 1.6 and rmax - 2.6 
(solid, black), r^nin = 2 and r^ax - 3 (dashed, red), r^in = 2.4 
and rmax = 3.4 (dotted, blue). Panel (d): the spectra from top to 
bottom are for /? = 2 (dotted, blue), /? = 1 (dashed, red), f3 - 
(solid, black), j6 = -0.7 (dot-dashed, magenta). 



i.e. larger h^/ rs, the hot spot under the source, where the variable 
part of the red wing is formed, is more extended and therefore 
the excess is less pronounced, see Fig.|5jb). 

In model S^^, the RMS spectrum is extremely sensitive to 
even small changes in r^in and r,nax, see Fig. EJc). For rmin > 
2.4, bending to the equatorial plane is too weak to give rise to 
a significant reflection component from larger r and the RMS is 
flat above 6 keV. 

Figure |5jd) illustrates changes of the RMS spectrum result- 
ing from the change of the radial profile of intrinsic luminosity, 
with the power-law modification parametrized by /3, defined in 
Sect. 12.21 Interestingly, less centrally concentrated profiles, with 
y8 > 0, preserve the shape of the RMS spectrum with only a 
slight change in the shape of the excess between 4 and 6; how- 
ever, the amplitude of the spectrum increases significantly with 
the increase of For /3 < 0, the RMS spectra flatten; further- 
more, for /3 < -0.5 their qualitative properties change, see the 
dot-dashed curve. In particular, no drop around 6.4 keV occurs, 
and the spectrum increases at £ > 10 keV. Note that such nega- 
tive values of y6 would characterise discs with non-zero stre ss at 
the marginally stable orbit, considered e.g. bv iKrolikl(ll999l) . 
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Fig. 6. Changes in the RMS spectrum corresponding to the in- 
crease of the number, nx, of the source positions contributing to 
the energy spectra. Panel (a): model S^^ with 5 = 0, rmin = 1.8 
and r,nax = 3. Panel (b): model C with = 4 and hmux - 20. 
Models in both panels assume a - 0.998 and Bobs - 35°. In 
both panels, the dashed (red online), solid (black) and dotted 
(blue) curves correspond to nx = 1,3 and 5, respectively. The 
dot-dashed (green) curve in panel (a) is for a single (nx = 1) 
primary source with the size Ar^ = 0.2. 



3.4. More complex arrangements of the primary source 

We consider here the RMS spectra constructed from the energy 
spectra, ft, formed by contributions from several positions of 
the source. A model with compact flares occurring simultane- 
ously at various random locations is an obvious example requir- 
ing such treatment. An equivalent effect should occur, regardless 
of the configuration of the emitting region, in any model involv- 
ing changes of location, for the increase of the size of time bins. 

The energy spectra in this (and the following) section are 
computed as a mixture of spectra from randomly generated po- 
sitions. In Fig.|6]the positions are drawn with uniform probabil- 
ity and in Fig.|2la) we illustrate the eff'ects of the change of the 
probability distribution. 

Figures |6ja)(b) show changes in the RMS spectra resulting 
from the increase of the number of source positions, nx, con- 
tributing to an individual spectrum, //. Obviously, deviation be- 
tween the average spectrum and each of the individual spectra 
decreases, with the increase of nx, which reduces both the RMS 
amplitude and the dependence of the RMS spectrum on energy. 

An essentially similar, but less pronounced, eff'ect occurs for 
an increase of the size, Ar^, of the primary source, see the dot- 
dashed curve in Fig.|6ja); emission of an extended source is ap- 
proximated by superposition of emissions from point sources lo- 
cated in the range [rs - 0.5 Arj, rs + 0.5Ars], so the energy spectra 
for Ars = 0.2 are formed by adding spectra from 3 adjacent po- 
sitions of the model S ^'^ grid. Again, deviations between spectra 
from more extended primary sources are smaller than in models 
with a single point source. However, diff'erences between spec- 
tra are more systematic, and larger, than between those resulting 
from nx > 1 with randomly generated rs and, therefore, the flat- 
tening effect is weaker 

Note that the change of the radial probability density (which 
effectively modulates the radial emissivity profile) has a diff'er- 
ent eff'ect on the RMS spectrum than the change of the radial lu- 
minosity profile. In particular, the dominating contribution from 
smaller rs (corresponding to 5 < 0) does not result in a qualita- 
tive change of the RMS spectrum (contrary to the modification 
with /? < 0), see Fig.lTja). 

Our model for the generation of active regions, with the 
same number of source locations contributing to each spectrum, 
and an implicit assumption of the same duration of emission at 
each position of the source, is certainly oversimplified; see e.g. 
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Fig. 7. Panels (a) and (b) are for model S^'^ with 0obs = 35°, 
a - 0.998, rmin = 1.8, r^ax = 3 and «x = 1- Panel (a): the spectra 
from top to bottom are for 6 = -1.5 (dotted, blue online), 6 = 
(dashed, red) and 5 = 1.5 (solid, black). Panel (b): The dotted 
(blue) curve is for constant intrinsic luminosity, L^^(rs), at given 
the solid (black) curve is for intrinsic luminosity generated 
from the Gaussian distribution with crL(rs) = 0.15L^^(rs); the 
dashed (red) curve is for intrinsic luminosity generated from a 
uniform distribution with s = 0.15. The inner panel (c) shows 
the RMS amplitude in a model involving changes of intrinsic 
luminosity as a function of the RMS variability without such 
changes; the dashed (red) and the solid (black) curves are for 
the uniform and Gaussian distribution, respectively, both with 
s = 0.15; see text for details. 



iPechacek et all (l2008l) for a more sophisticated modeling of the 
generation of flare-like features in random processes. Our main 
purpose here is to illustrate flattening of the RMS spectrum, 
whi ch is a generic trend corresponding to the increase of rix- See 
also ICzernv et al.l ([20041) for the RMS spectra from their model 
of a spotted disc, which are much flatter than those made here 
with small rix (however, their model neglects effects of transfer 
from the X-ray source to the disc, which is a qualitative differ- 
ence to our model S '^^). 

3.5. Changes of intrinsic luminosity 

Finally, we take into account changes (in time) of intrinsic lu- 
minosity. The ability to describe the entire variability pattern in 
terms of relativistic effects is an attractive property of GR mod- 
els, however, assumption of a strictly constant luminosity seems 
rather unrealistic. Results of the previous sections, though, could 
be applied directly to modeling a system with the time-scale of 
intrinsic variations much shorter than the time-bin length (so 
that each bin probes emission averaged over various luminosity 
states). 

Figure |3b) shows the RMS spectra for model 5''^^, combin- 
ing changes of the radial location with changes of the intrinsic 
luminosity (at each given r^). The location is randomly gener- 
ated (with (5 = 0) and for each we generate a random value of 
luminosity. Thus, the generated sequence of basic energy spec- 
tra typically contains several contributions from a given location, 
each with a different normalization. 

It is interesting to notice that the results depend on the as- 
sumed distribution of luminosities. In order to illustrate this, we 
construct the RMS spectrum using: 

(i) a Gaussian distribution, P(L) oc exp[-0.5(L-Lo)^/o"L]' where 
Lo = LPT(rs) and (Tl = sl'^ir,); 

(ii) a uniform distribution in the range [Lo(l - -^), Lo(l + s)], with 



In both cases, 5 is a constant number (independent of rj, thus, 
both the deviation for distribution (i) and the amplitude of vari- 
ations for distribution (ii) change slightly with r^. In general, 
changes of intrinsic luminosity increase the amplitude of the 
RMS spectrum, however, the Gaussian distribution leads also to 
flattening, while the uniform distribution does not change the 
shape of this spectrum with respect to that from the model with 
constant luminosity. 

The relevant property underlying the above difference is il- 
lustrated in Fig.lTJc). The figure shows the relation between the 
RMS computed with (RMSgr+l) and without (RMSqr) changes 
of intrinsic luminosity for a simple model using two random 
quantities, Xf and xl. xf is generated using a uniform distribu- 
tion and represents the received flux in GR models with con- 
stant luminosity, xl represents the intrinsic luminosity and is 
generated using a Gaussian (the solid curve) or uniform (the 
dashed curve) distribution, as defined above, both with s = 
0.15. RMSgr is constructed using the generated values of Xf, 
the increase of the amplitude of variations of Xf yields larger 
RMSgr. RMSgr+l is constructed using the products of Xf and 
XL. Note that the uniform distribution of xi leads to a linear 
relation between RMSgr+l and RMSgr, which is reflected in 
the unchanged shape of the spectrum in Fig. |2lb), while for the 
Gaussian distribution this relation flattens at small RMSgr. 



4. Application to l\/ICG-6-30-1 5 

4.1. Tlie Suzaku observations 

Because of the wide bandpass coverage of energies provided by 
detectors on-board Suzaku, it is currently the most suitable X- 
ray observatory for testing relativistic reflection models. In par- 
ticular, crucial information comes from the Hard X-ray Detector 
data, allowing us to measure the Compton reflection hump si- 
multaneously with the Fe Ka line. MCG-6-30-15 was observed 
three times by Suzaku in 2006 January in a state typical for 
this source, regarding the mean X-ray flux and variations of the 
flux by a factor of ~ 2 in a few ksec (see Miniutti et al. 200?!). 
Crucially for our analysis, Miniutti et al.l 12007 ) find for these 
observations that (1) the average profile of the Fe Ka line has 
a relativistic shape, with the pronounced blue peak and the red 
wing extending down to ^ 4 keV, consistent with that from the 
long XMM observation in 2001 (Fabian et al. 2002); (2) the av- 
erage relative normalization of the Compton hump, given by the 
usual parameter R {R = 1 corresponds to the strength of reflec- 
tion expected from a reflector subtending 27r sr at an isotropic 
primary source), isR ^ 3.8, consistent with previous estimations 
from Beppo-SAX (Ballantyne et al. 2003), and rises io R ^ 5 
in the low flux spectra; (3) on a time-scale of tens of ksec, 
spectral changes are consistent with the phenomenological, two- 
component model. 

We use the data from three observations of MCG-6-30-15 by 
the Suzaku satellite on 2006 January 9-14 (143 ksec exposure), 
23-26 (99 ksec) and 27-30 (97 ksec). The total on-orbit time is 
about 780 ksec. For data reduction, we have used the HEADAS 
6.5 software package provided by NASA/GSFC. We determined 
the RMS amplitude, and its error, as a function of energy, in 
the standard manner (e.g. lVaughan et aLll2003b; the definition is 
analogous to equation (|3]l but contains an additional term related 
to the count rate error). In this paper we use the RMS spectra for 
two time-bin sizes, approximately 16 ksec and 131 ksec (more 
specifically 2'^* sec and 2'^ sec), with approximately 7 and 50 
ksec of exposure time per bin, respectively. The RMS spectra 
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Fig. 8. All curves show the RMS spectra for a model with a con- 
stant Compton reflection (derived with pexrav), a constant line 
at 6.4 keV with EW - 30 eV (and with cr = 0.1 keV, match- 
ing the width of the data energy bin) and a varying power-law 
component. Panel a: the RMS spectra are computed with var- 
ious amplitudes of the power-law variability, but all have the 
same R = 0.15, approximately consistent with the EW of the 
line. Panel b: the simulated RMS spectra, with R = 4, com- 
pared with the Suzaku data. The upper (red) and lower (magenta) 
points show the RMS spectra for the Suzaku observations with a 
time bin size of 16 ksec and 131 ksec, respectively. 

for 16 ksec and 131 ksec are shown in Fig.[8lb) and all panels of 
Fig.|9]by the upper and lower set of points, respectively. 

The RMS spectra indicate strong decrease of variability at 
< 1 keV, however, modeling of this band is beyond the scope of 
this paper The decrease most likely results from the presence of 
a soft X-ray excess, which is clearly seen in the average spectra 
of MCG-6-30-15 as well as in the spectrum of the non- varying 
spect ral component (see e.g. figs. 3 and 22 in lVaughan & FabianI 
12004 . The origin of such soft excesses, typically observe d 
in AGNs, is not clear (see e.g. ISobolewska & Don3 l2007h . 
Remarkably, one of the possible explanations, involving the rel- 
ativistically blurred photoionized disc reflection (where the soft 
excess is expl ained as being composed of many broad lines, see 
ICrummv et al..,2006i) . is consistent with the scenario discussed 
in Sect. 1531 

4.2. Non-relativistic reflection 

Before applying the GR models, we briefly comment on the con- 
tribution of reflection from distant matter, which for some other 
objects is considered as the explanation for RMS spectra similar 
to these derived in MCG-6-30-15 (e.g. Terashima et al. 2009). 
In MCG-6-30-15, however, such a component should have a 
marginal effect, as indicated by the very small contribution of a 



narro w Fe Ka line at 6.4 keV (e.g. llwasawa et al.ll996l : lLee et al.l 
l2002h . 

For the Suzaku observation, the equivalent width (EW) of 
the n arrow 6.4 keV line is approximately 30 eV jMiniutti et al.l 
120071) . suggesting that the distant reflector subtends a small solid 
angle at the central X-ray source with R - 0.15 being a likely 
value characterising the strength of the accompanying reflection 
hump. In Fig. |8^, we show the RMS spectra for a simple non- 
relativistic model involving a variable power-law and a constant 
reflection component with R - 0.15 and a narrow Fe Ka line 
with EW=30 eV (both EW and R are determined with respect to 
the average power-law component). The RMS spectra are com- 
puted for various levels of variation of the power-law contin- 
uum; the index, F, of the power-law remains constant while its 
intensity changes. The reflection spectrum is derived using the 
pexrav model, with the inclination angle of the reflector fixed 
at 30°. The constant line at 6.4 keV is computed with a rather 
large width, cr = 0.1 keV, for a clear comparison with the ob- 
served RMS spectra (which are derived with the energy bins 
AE = 0.2 keV around 6.5 keV). The constant line at 6.4 keV, 
with EW = 30 eV and cr = 0.1 keV, is also included in all mod- 
els in Figs.[8l:b),|9]and[l0] 

For a constant spectral component, with the fixed contribu- 
tion to the average spectrum, the strength of the related signa- 
ture in the RMS spectrum depends on the overall level of vari- 
ability. For RMS = 0.3, both a constant hump with R - 0.15 
and a line with EW - 30 eV would lead to rather pronounced 
declines in the RMS spectrum (see the top curve in Fig. [8^). 
However, if another spectral component dilutes the variability, 
as should be the case in MCG-6-30-15, the strength of these de- 
clines decreases. Specifically, RMS = 0.05-0.15, derived from 
the data above 10 keV, makes the hump with R - 0.15 negligi- 
ble. For RMS = 0.15-0.2, as observed around 6 keV, the line 
with EW = 30 eV has a noticeable, but rather minor, eff'ect. 

The observed shape of the RMS spectrum can be approx- 
imately reproduced, above 3 keV, in a simple non-relativistic 
model, involving a constant reflection component with R 3 
(see Fig.lSJ)). Such an apparently unphysical (for reflection from 
a distant matter) value of R could result from the effect of the 
light travel time in a source observed in low luminosity state (as 
for the Suzaku observation of NGC 405 1 ; Terashima et al. 2009) 
or from GR effects affecting the primary emission. In the latter 
case, more relevant for MCG-6-30-15, the mechanism is simi- 
lar to model S^^. Namely, if generation of X-rays occurs within 
a few Rg from a rapidly rotating black hole, a face-on observer 
would directly receive reduced primary radiation, which would 
be strongly focused along the equatorial direction to the distant 
material, giving rise to a strong non-relativistic reflection com- 
ponent. In MCG-6-30-15 this scenario seems to be ruled out by 
the small EW of the narrow Ka line. Moreover, modeling the 
RMS spectra for various time-bin sizes seems to require slightly 
different values of R (approximately 3 and 4 for 16 ksec and 131 
ksec, respectively), which is inconsistent with a simple reflec- 
tion model (where R should not depend on the time-bin size). As 
an illustration, the RMS spectra in Fig. |8lb) are computed with 
R - 4. While the lower curve roughly matches the spectrum for 
131 ksec, the upper curve strongly underpredicts the variability 
above 10 keV, as compared to the spectrum for 16 ksec. 

4.3. GR models 

The RMS spectra predicted by the GR models for some ranges 
of parameters are qualitatively inconsistent with those observed, 
as discussed in Sect. [3] Thus, the observed data completely rule 
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Fig. 9. Fits of the Suzaku RMS spectra with the GR models. In all panels, the upper (red online) and lower (magenta) points show 
the Suzaku RMS spectra for 16 ksec and 131 ksec, respectively. All models involve a constant line at 6.4 keV with the EW = 30 
eV and cr - 0.1 keV; all models assume a - 0.998 and flobs - 35°. Panel (a): models involving vertical motion of the primary 
source. The dashed (green online; shifted vertically by -0.4) curve is for model A with h^^^ - 2.2 and /imax - 10. The solid (blue; 
shifted vertically by -0.07) and dotted (black; shifted vertically by 0.01) curves are for model C, with /jmin = 4 and /zmax - 20, for 
nx = 1 and 3, respectively. Panel (b): the solid (blue) and dotted (black) curves are for nx = 3 and 5, respectively, in model S^'^ 
with Ttnin - 1.6, Tmax = 3, 5 = - 1 .5, (Tl = 0.2L^^ (and j6 - 0). The dot-dashed (black) curve is for the same parameters as the solid 
curve, except for (larger) rmax - 3.4. The dashed (green) curve is for r^in = 1.4 and rmax = 3; nx = 3, j6 = 0, 5 = 0, CTg - 0. Panel 
(c): model 5 with r^^^ = 1.8 and r^ax - 3, /3 = 1 (5 = and ctl = 0); nx - 2 and 3 for the solid (blue) and dotted (black) curve, 
respectively. 



out model A (for completeness of our discussion, the best fit with 
this model is presented in Sect. 14.3.11 ). model S with > 4 and 
model C with < 2, cf. Figs. [1] and [3] 

Some details of the RMS spectra in the GR models are sen- 
sitive to the value of Oobs- In our analysis we consider the in- 
clination angles between 0obs - 25° and 45°. We recall that 
the disc-line fits to the ASCA observations of MCG-6-30-15, 
in which the line seems to extend only up to ^ 6.8 keV, yield 
Oohs - 30°. However, higher quality data indicate that the line 
extends beyond 7 keV (but is in dented by two absor ption edges) 
and recent fits to bo th the XMM dLarsson et alj2007l) and Suzaku 
jMiniutti et alj|2007ft observations indicate a larger value, flobs - 
40°. 

In models A and C we adjusted the amplitude of the RMS 
spectra to match the observed data, which is the usual procedure 
in spectral modeling. In model S^^ we did not treat the ampli- 
tude as a free parameter. Instead, we attempted to reproduce self- 
consistently both the shape and the amplitude of the RMS spec- 
trum by modifying the auxiliary parameters of the model, i.e., 
(Tl and 6 or /3. 

The model spectra best fitting the data are shown in Fig. |9l 
Our formally best fitting case, shown by the solid curve in Fig. 
Igtb) has the (reduced) xl ~ 2 (with 24 d.o.f.) for £ > 3 keV; 
the remaining fits in Fig. |9] yield > 3. Clearly, the fits are not 
acceptable statistically. On the other hand, the basic trends are 
reproduced by model S ^'^ reasonably well and we suspect that 
the quality of the fits could be improved when additional phys- 
ical effects, neglected in this paper (e.g. ionization), are taken 
into account. 

A further indication of the incompleteness of our models is 
given by the systematic discrepancies below 4 keV. Between 2 
and 4 keV, the RMS spectra predicted by the GR models are 
typically more convex than these observed, similarly to the case 



of non-relativistic models (cf. Fig. [8j)). Only model S^^ with 
r,nin < 2 (producing a strongly redshifted reflection component), 
yields the RMS spectra with approximate agreement down to 2 
keV. 

Studies of strong gravity effects in the time-averaged spec- 
trum of MCG-6-30-15 usually concentrate on the data above 
3 keV, as the spectrum at lower energies is strongly affected 
by a warm absorber (while at larger energies the absorption ef- 
fects are considered to be unimportant; see e.g. iMiniutti et al.l 
12007'). We also do not attempt to model the RMS spectra below 
3 keV, as changes of the absorber may determine variability in 
this band. On the other hand, some controversies remain as to 
whether the absorber does exhibit strong changes, see Sect. 15.61 
and we point out that a weakly varying absorption should not 
affect the RMS variability significantly. Considering the lack of 
absorber variability suggested by some previous studies, we note 
that the deficiency of GR models below 3 keV results from the 
small contribution to the total spectrum of the radiation reflected 
from neutral matter at these energies. Then, ionization of the in- 
nermost parts of the disc (which is most likely in model 5^^^, 
see Sect. 15.51 ) may be relevant to extend these models to lower 
energies. 

4.3.1. Models A, C 

In models A and C we considered < 20. We were not able 
to find the RMS spectra matching both the shape and the nor- 
malization of those observed. Then, we focused on reproducing 
the shape of the RMS spectrum and we adjusted the amplitude, 
which is treated as a free parameter only in this section. We pre- 
sume that moderate changes of the RMS normalization could re- 
sult from a departure from the (vertically constant) distribution 



A. Niedzwiecki & T. Miyakawa: GR models of the X-ray spectral variability of MCG-6-30-15 



11 



of the intrinsic luminosity assumed in these models. In general, 
the model RMS spectra with a (qualitatively) consistent shape 
have larger normalization than those observed, therefore, we ne- 
glected variations of intrinsic luminosity in time, which would 
lead to a further increase of the RMS amplitude. 

The dashed curve in Fig. |9la) shows the RMS spectrum in 
model A, with Sobs - 35°, /;,nm = 2.2 and h,^^^ = 10, best match- 
ing the observed shape; the model spectrum is shifted down by 
ARMS = -0.4. An excess variability in the energy range of the 
Fe Ka line, robustly predicted by model A, is clearly inconsis- 
tent with the observed data. We do not expect that modifications 
of our model, specifically those related to ionization of the disc, 
could reduce this discrepancy (see Sect. 15.5b and thus we rule 
out this model. We also note that the observed level of the RMS 
variability can be achieved in model A with h^i„ > 3, but in such 
a case the model RMS spectrum is much flatter than observed. 
On the other hand, the stronger contribution from small h^, e.g. 
for /imin < 2.2, yields a larger excess around 6.4 keV. 

Model C predicts, for Oobs > 30°, a sharp drop at the max- 
imum of the blue peak, cf. Fig. 1(c), which cannot account for 
the observed, broader depression in 6-7 keV range. However, 
for flobs = 35°, the drop occurs around 6.8 keV and including an 
additional narrow line at 6.4 keV we achieve an overall form of 
the depression between 6 and 7 keV which can imitate the ob- 
served one. Then, we found the set of parameters (Oobs = 35°, 
hmin = 4, /imax = 20; the solid curve in Fig.|9^) for which devia- 
tions between the predicted and observed shapes are small. The 
predicted RMS amplitude only slightly exceeds the observed 
value, the model spectrum is shifted down by ARMS = -0.07. 
The dotted curve in Fig. |9la) illustrates the extrapolation of the 
model for larger time-bin sizes, obtained with «x > 1 ■ 

Model C requires a fine-tuning of parameters to approxi- 
mately explain the observed fractional variability. For the in- 
clination angles higher or lower than 0obs - 35°, discrepancies 
around the energy of the Fe Ka line are more significant. For the 
range of heights different to ^ 4-20, the RMS spectrum be- 
tween 3 and 20 keV is too flat. See Sect. |5.3| for a further, critical 
discussion of the fits with model C. 

4.3.2. Model 5 

Both the shape and amplitude of the observed RMS spectrum 
can be approximately reproduced with a > 0.98 and the distance 
of the primary source varying in the range extending down to 
at least 2. For smaller spins, or larger distances, the model 
RMS spectrum is flatter than observed. Models with nx > 1 
strongly favour the maximum value of a = 0.998. The major 
challenge in modeling the observed RMS spectrum results from 
the excess related to changes of the red wing, which is not seen 
in the data. A strong ionization of the disc, at the site of the 
formation of the variable part of the red wing, may be relevant 
in reducing this discrepancy, see Sect. 15.41 In the simplest case 
(with, in particular, a neutral and untruncated disc) considered 
in this paper, fits with model S^^ favour r^in < 2, for which 

(1) the excess is more extended and less pronounced, making a 
relatively small deviation between the model and the data, and 

(2) the significant contribution from the smallest allows us to 
reproduce the RMS spectra down to 2 keV. 

We attempted to reproduce the RMS spectra for both 16 and 
131 ksec with the same set of parameters, differing only by the 
value of «x. In our procedure, we first adjusted r^i^, '"max, "x, <5, 
(Tl and flobs to fit the 16 ksec spectrum, then we computed the 
spectra for larger nx and compared them to 1 3 1 ksec. The differ- 
ences between the amplitudes of the RMS for consecutive nx are 
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Fig. 10. RMS spectrum for model S^^ with an excess around 
5.5 keV, as revealed by some previous observations (see text). 
Parameters = 1 .2, 5 = 0, ctl = and nx - I are adjusted to 
match the Suzaku RMS spectrum (for 16 ksec) at £ < 5 keV and 
E>6 keV. 

rather large, which is most challenging in attempts to reproduce 
both time-bins simultaneously (as noted before, the underlying 
assumption of identical nx for each time bin is most likely over- 
simplified). 

Our best fits, with a = 0.998, 0ob,s = 35°, r^u, = 1 .6, r^ax = 3, 
ctl = 0.2L^^, 6 = -1.5, nx = 3 and 5 for 16 and 131 ksec, re- 
spectively, are shown in Fig.|9jb). Disregarding the excess vari- 
ability of the red wing (with the plausible impact of ionization in 
mind) we can reproduce the observed RMS spectra for a broader 
range of parameters. Example fits, achieved using p instead of 
(Tl and 6, with Oobs - 35°, r^in - 1-8 and r^ax = 3,yS = 1, nx = 2 
(for 16 ksec) and 3 (for 131 ksec) are shown in Fig.|9lc). 

The statistical quality of the fits is poor, as noted above, 
therefore, the usual criterion for constraining the confidence lim- 
its in the model parameters, using the increase of;^'^, seems in- 
adequate. Qualitatively, r^in - 2.2 seems to be the upper limit 
for Tniin; for larger rmin, the RMS spectrum is too flat, see the 
bottom curve in Fig. 5(c). Similarly, we can reject submaximal 
(< 0.98) values of a, cf. Fig. 4(a). The model with r,^in -1.6 and 
'"max = 3 is slightly favored in our fits; Fig. |9jb) shows trends 
corresponding to the extention of this range. In general, mod- 
els with r^ax > 3 give a more significant excess with respect 
to the data below 5 keV, see the dot-dashed curve. In turn, for 
our best fit with r^i„ = 1.4, variability above 10 keV is slightly 
underpredicted; notably, this fit (the dashed curve) involves no 
modification of our basic S^'^ model, i.e. ctl = 0, 6 - and 
yS = 0, so the intrinsic luminosity follows exactly LF^. See Sect. 
15. 41 for further discussion of r^i„ and r^ax- 

4.3.3. Excess variability between 5 and 6 keV 

As we discuss in Sect. 15.51 the excess variability of the red 
wing, which is the major prediction of model S^'^ ruled out by 
the Suzaku data, is likely to be reduced by ionization of the 
disc under the source. Interestingly, however, the RMS spec- 
tra from some previous observations of MCG-6-30-15 (see 
Matsumoto et al. 2003; Ponti et al. 2004) did reveal a strong en- 
hancement of the RMS variability between 5 and 6 keV. Apart 
from this band, those RMS spectra are consistent with other ob- 
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servations of MCG-6-30-15. Presumably, the geometry of the 
X-ray emitting region does not change significantly between the 
various observations and the presence or lack of an excess is re- 
lated to changing ionization state of the innermost disc. 

Then, we examined conditions in which model S produces 
an excess between 5 and 6 keV but matches the Suzaku RMS 
spectra at other energies. An example of such an RMS spectrum 
for 0obs - 40°, r^in = 2.2 and r^ax = 3.2 is shown in Fig. [TOl 
We find that, in general, a relatively high inclination angle, ^ 
40° - 45°, is required to explain the excess between 5 and 6 keV 
by a varying redshift of 6.4 keV photons - for smaller 0obs the 
excess occurs at lower energies. 

The amplitude of the e xcess in Fig . lTOlis a pproximately con- 
sistent with that found by iPonti et al] ( 120041) in the XMM data. 
However, we were not able to reproduce a stronger excess, with 
RMS ~ 0.4, as reported by Matsumoto et al. (2003) for the 
ASCA observation. Such a strong excess can be explained by 
an enhanced contribution of Fe photons from the spot under the 
source, which could occur if ionization of that area of the disc is 
taken into account, namely, it could be attributed to intermediate 
stages of ionization (where Auger destruction cannot operate) 
giving rise to an intensity of the Fe photons much higher than 
from a neutral medium (e.g. IZvcki & Czernvlll994l) . For these 
ionization states, the rest energy of Fe photons is around 6.9 keV, 
and the range of inclination angle relevant to reproduce the ex- 
cess would be 6»obs - 35° - 40°. 

Interestingly, during the XMM observation analysed by Ponti 
et al. (2004), the source luminosity was ~ 2 times lower than the 
typical luminosity observed in MCG-6-30-15. This can support 
the above scenario, with the change from complete to interme- 
diate ionization of a hot spot. In case of two ASCA observations 
used by Matsumoto et al. (2003) to reveal enhanced variability 
between 5 and 6 ke V, the average luminosity was not much lower 
than typical, however, both observations contain prolonged deep 
minimum states, with the luminosity dropping to the lowest level 
observed in that object. 

We point out that an excess between 5 and 6 keV cannot be 
produced by simple geometric effects in models A or C; actually, 
the latter predicts an opposite property, i.e. reduced variability 
around 5.5 keV. 



(2) Relativistic distortions of the Fe Ka line give rise to excess 
variability which is the major challenge in applications of GR 
models to observed data. In models involving vertical motion of 
the X-ray source, the excess variability, between 6 and 7 keV, is 
related to changes of the blue peak. In models involving changes 
of the radial distance, low above the disc surface, the excess is 
related to changes of the red wing and occurs between 4 and 6 
keV. In the latter class of models, however, ionization of the hot 
spot under the source is very likely and should reduce the excess 
variability, see Sect. 15.51 

(3) The results of this paper illustrate the maximum strength of 
signals in the RMS spectra caused by the GR effects. Such strong 
signals could be observationally seen only in the RMS spectra 
determined with time-bin sizes not exceeding the time scale of 
the change of position of the X-ray source. For longer time bins, 
signatures of GR effects are weaker (i.e., the RMS spectra are 
flatter); some implications of this are discussed in Sect. 15.21 

(4) Variability observed in MCG-6-30-15 is inconsistent with 
models assuming vertical motion of the source. The model with 
a source on the symmetry axis is ruled out and the model by 
MF04 involving a specific pattern of rotation around the sym- 
metry axis is disfavored as discussed in Sect. 15.31 

(5) The model wit h chang es of the radial distance in the Kerr 
metric, proposed in iNZOSl offers the most likely explanation of 
the observed properties (among the GR models); See Sect. 15.41 
for further discussion of this model. 

(6) Some details of the RMS spectra are sensitive to 9obs, and 
they may provide constraints on inclination, independent of the 
modeling of the line profile. Interestingly, our GR fits to the 
RMS spectra from Suzaku observations of MCG-6-30-15 fa- 
vor 0obs = 35° - 40°, consistent with the disc-line fits to the 
Suzaku data. The light-bending (MF04) model requires exactly 
6»obs = 35°. The best fit with model 5 also has Ooh^ = 35°, how- 
ever, this model allows for a broader range of Oobs, especially 
when excess due to variability of the red wing is disregarded. 
The range of 0obs = 35° - 40° would also be relevant in repro- 
ducing the excess variability between 5 and 6 keV, which was 
reported for some previous ASCA and XMM observations, in a 
model involving a highly ionized hot spot under the hard X-ray 
source. 



5. Summary and discussion 

5.1. Conclusions 

This paper focused on modeling the RMS spectra, which are 
commonly used to estimate the fractional variability amplitude 
in AGNs. We thoroughly studied the RMS spectra predicted 
by models relating spectral changes to varying amounts of GR 
eflBCts, which distort the X-ray radiation generated in central 
parts of black-hole accretion flow. We applied our GR models 
to MCG-6-30-15 where the bulk of the reprocessed radiation is 
supposed to come from the inner accretion disc. Our main con- 
clusions are summarized as follows: 

(1) The Compton reflection hump, reprocessed in the accretion 
disc from the continuum emission generated at a small, varying 
distance from a black hole, is typically less variable than the pri- 
mary continuum. This leads to a decrease of the RMS spectrum 
with increasing energy, above 3 keV (for neutral disc material). 
A remarkable exception involves a source corotating with the 
disc and changing height (as in the model of MF04) not far above 
the disc surface - in this case the fractional variability increases 
with energy. 



5.2. Variability time-scales 

The energy spectra determined for increasing lengths of time- 
bins should consist of contributions from an increasing number 
of locations of the primary source. Then, the RMS spectra for 
longer time bins should be flatter, see Sect. 13. 41 Below we briefly 
discuss the resulting implications for time-scales of various pro- 
cesses involved in generating the observed variability. 

The reduced variations of the reprocessed component are 
typically assessed over ^ 10 ksec. We emphasize that the actual 
time-scale of the change of the position of the hard X-ray source 
should be at least of this order of magnitude, otherwise the RMS 
spectra would be much flatter than observed. Then, significant 
changes of the X-ray flux on a time-scale of ~ 1 ksec which 
are observed in MCG-6-30-15 should result from intrinsic varia- 
tions (with small changes of geometry) of the hard X-ray source. 
The point-to-point RMS spectrum with 1 ksec time bins which 
measures such short time-scale varia bility, found to be much flat - 
ter than the standard RMS spectra in IVaughan & FabianI (|2004|) . 
supports the scenario with changes of intrinsic luminosity domi- 
nating the short time-scale variability (by definition of the RMS 
spectrum, changes of luminosity, without the change of spectral 
shape, produce a flat RMS spectrum). 
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In Sect. 14.3.21 we note that the RMS spectra for both 16 ksec 
and 131 ksec can be reproduced with the same set of parameters, 
in a model involving a larger number of the source positions for 
the larger bin size. However, the number of positions in our fits 
for 131 ksec is only ~ 2 times larger than for 16 ksec, while a 
factor of 8 would be expected for a simple scaling of the number 
of positions with the bin length. Such a simple, linear scaling 
of the number of source locations is ruled out, as ^ 10 random 
locations, implied for 131 ksec, would yield an almost flat RMS 
spectrum. Then, the change of geometry on a time scale of ~ 100 
ksec should obey some systematic pattern, involving at most a 
few localised regions dominating the total emission, with the du- 
ration of several tens of ksec. Emission averaged over the whole 
range of locations of the source is probably probed on a much 
longer time-scale, > 1000 ksec, e.g. by the RMS spectra derive d 
from long-term RXTE monitoring, see iMarkowitz et al.l (l2003l) . 
which are indeed flatter than these obtained with 10 - 100 ksec 
bins. 

Finally, note that the time-scales for the change of the source 
location, ~ 10 - 100 ksec, required by the above arguments are 
orders of magnitude longer than the dynamical time-scale for 
the innermost region, ^ 100 sec, which is rather puzzling. On 
the other hand, the power-spectra indicate that most of the vari- 
ability in MCG-6-30-15 indeed occurs on time scales > 10 ksec 
and a large-magnitude varia bihty may extend even to time-scales 
longe r than 100 ksec (see IVaughan et alj 1200^ lUttlev et all 
I2OO2I) . Then, in the context of GR models, most of the power 
contained in the X-ray light curve of MCG-6-30-15 would re- 
sult from changes of the geometry of the X-ray source. In turn, 
variability dominated by changes of intrinsic luminosity should 
correspond to frequencies higher than the break frequency ( 10 
Hz, above which the magnitude of variability drops). 

5.3. Models involving vertical motion 

Focusing of radiation produced close to the symmetry axis, to- 
ward the disc, cannot explain the detailed properties of MCG- 
6-30-15, despite qualitative arguments for the relevance of this 
mechanism. Significant changes of the blue peak corresponding 
to changes of the source height are the major cause of this model 
failure. 

This discrepancy can be reduced in a modification of the 
model by MF04, where the observed RMS spectrum can be ap- 
proximately reproduced with the height changing between ^ 4 
and 20Rg. We emphasize again that predictions of this light- 
bending model are related primarily to its specific kinematic 
assumptions rather than the light-bending itself. We disfavor 
this model on the grounds that it requires large elevations of 
the source above the disc, where the kinematic assumptions are 
rather arbitrary. At lower elevations (where the crucial assump- 
tion of this model - i.e. rotation of the source with Qk(p,s), see 
Sect. 12.11 - has some physical motivation), the model predicts a 
variability pattern that is ruled out by observations. Even with 
the required range of h^, the model is rather unlikely as a de- 
cline of the fractional variability around 5.5 keV, systematically 
predicted by the model, has not been confirmed by any observa- 
tion of MCG-6-30-15 (actually an opposite property has been 
reported, see Sect. |4.3.3t . Another major shortcoming of this 
model was pointed out by iNZOSi For the range of heights rel- 
evant to reproduce the observed variability (> 4-Rg), the model 
cannot explain the pronounced red wing of the relativistic Fe 
line, the feature that originally motivated the development of this 
model. 



5.4. Model with bending to the equatorial plane 

We strongly favor the model with radial motion in the Kerr met- 
ric for a description of the central region of MCG-6-30-15. Note 
that both of the main properties required by this model, i.e. rapid 
rotation of the black hole and very small distance of the X-ray 
source, have some independent support in the observed, time- 
averaged Fe line profile. Our condition for the value of the spin, 
a > 0.98, is consistent with the formal c onstraint on the spin 
parameter, a > 0.987 jBrenneman & Reynolds 2007), derived 
from the line profile under the assumption that the Fe photons 
are emitted only at distances larg er than r^^ia). Regarding the 
second condition (small distance), iNZOSi point out that location 
of the X-ray source within a few central Rg is needed to explain 
the pronounced red wing of the line profile; the irradiation of the 
inner disc by a source located at further distances is too weak 
and the red wing is much weaker relative to the blue peak than 
in the observed profile. 

In the relevant range of distances, the directly observed flux 
of primary X-rays is extremely sensitive to the location of the X- 
ray source, while the reflected radiation changes a little, and the 
observed variability effects can be reproduced by rather small 
changes in the geometry of the X-ray emitting region. The range 
of distances between approximately 1.6/?c, and 3Rg is slightly 
preferred in our applications of this model to the Suzaku data. 
However, we emphasize again that some features of our model 
are oversimplified, most importantly our method for the genera- 
tion of active regions and their luminosities as well as our neglect 
of ionization, and more realistic modeling could yield a more ex- 
tended range of locations. 

Remarkably, for a - 0.998, the max imum of the dissipation 
rate given by the Page & Thornd (11974 ) formula occurs exactly 
at the lower end (i.e. l.bRg) of this range. Furthermore, accord- 
ing to that formula with a - 0.998, 35 per cent of the total power 
dissipated in a Keplerian disc is released within the central 3Ra. 
Then, our condition for the location of the hard X-ray source, 
although seemingly rather extreme, is not implausible. 

An independent constraint on the range of distances comes 
from the normalization of the reflection component in the time- 
averaged spectrum. We have computed the reflection parame- 
ter, R, in our models as the ratio of the total fluxes in the re- 
flected component in the GR model and in the pexrav model 
with the same flux of the primary power-law. Our fits with model 
S^^, shown in Fig. |9] have R ^ 4, within the con fidence limits, 
R - 3 .8 + 0.7, found for Suzaku observations in iMiniutti et al.l 
(l2007h . A stronger contribution of primary hard X-rays from 
smaller radii (e.g. for r^i„ < 1.6, r^ax = 3 and S < 0) would 
yield, however, a larger value of R, inconsistent with the data. 

In our investigation of the model we focused on the case 
with /is/^s < 0.1. For /Js/^s » 1, i-e. in the lamp-post regime, 
the bending to the center dominates over the bending to the 
equatorial plane, leading to qualitatively different effects, as de- 
scribed above. We did not investigate the case of intermediate 
/zs/^s, as results in this regime strongly depend on (very uncer- 
tain) kinematic assumptions. Note that even for relatively small 
/js/^s = 0.2, the difference between the azimuthal velocities re- 
sulting from the different assumptions discussed in Sect. |2.1| ex- 
ceeds 0.1c. 

Note that occasionally an additional source of primary X-ray 
emission may be present in MCG-6-30-15, which would lead to 
a more complex temporal behaviour than estimated in this pa- 
per. Flaring releases of a part of the X-ray emission may occur at 
rather large heights (> 100/?g, note that GR effects are negligible 
at these distances). We note the following hints for such episodic 
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releases. MCG-6-30-15 typically shows no evidence that varia- 
tions in the Fe line track the continuum variations, which sup- 
ports the scenario with X-rays released and reprocessed in the 
very central region. However, during very strong flares (in the 
light-curve), the line appears to respond to changes of the X- 
ray flux and lags of a few ksec are observed between the flare 
and the a ppearance of an enhanced Fe line ( Ponti et aTl 120041: 
iNegoro et al. 2000). A likely explanation is given by the strong 
X-ray flare at a height of several hundred Rg. In agreement with 
such scenario, the (delayed) Fe line became narrower in time 
as if it was coming from more distant regions of the disc, in a 
manner consistent with a response to the flare that occurred at a 
height of several hundred Rg. 

Finally, we note a caveat against the specific scenario consid- 
ered in this paper, with hard X-rays produced by (at most) sev- 
eral compact sources corotating with the disc. A quasi-periodic 
modulation of the primary emission should appear in such a case 
on the Keplerian time scale, ~ 100 sec, however, such a sig- 
nal is not observed in MCG-6-30-15 (see lZvcki & Niedzwieckil 
Hoollfor details). This may imply that the hard X-ray source has 
a more continuous spacial distribution, e.g. forms an extended 
corona covering a part of the disc surface or a small hot torus 
replacing the innermost disc. If such an extended, hot plasma is 
located in a very small central region and exhibits small changes 
of geometry, the mechanism described in Sect. l3.3l should lead to 
properties qualitatively consistent with t hese d erived for model 
5 - see iNiedzwiecki & Frankiewic3 (l2007b for preliminary 
results of the model with a small, shrinking/expanding inner 
torus. 

5.5. Ionization 

The neglect of the ionization of the disc surface is the major 
shortcoming of our model. Determination of the ionization struc- 
ture requ ires rather complicate d radiative transfer computations 
(see, e.g.. lR6zariska et alj2002l) which are currently not included 
in our model. For MCG-6-30-15, it seems to be well established 
that the blue peak of the Fe Ko- line arises from neutral Fe at a 
distance of ^ lORg (e.g. Vau ghan & Fabian 2004). However, the 
inner parts of the disc, where the red win g is formed, may be 
highly ionized (cf. iBallantvne et al.ll2003h . The irradiation pat- 
tern in model S should yield a non-uniform ionization struc- 
ture of the disc, with a strongly ionized inner region and neutral 
outer parts, as for the relevant range of distances, Ts = 2 - 3, the 
flux irradiating the spot immediately under the source is ~ 4 or- 
ders of magnitude larger than that irradiating the disc at ~ IQRg. 

A precise estimation of the impact of ionization effects on 
the Fe line variability is difficult without detailed computations, 
as a compact X-ray source above the inner disc should produce 
a range of ionization states over the surface of the disc and the 
profile of ionization should change in response to the change of 
location of the source. Qualitatively, however, we may expect 
that fits of the observed RMS spectra would be improved with 
these effects taken into account, in particular the major discrep- 
ancy between model S^'^ and the Suzaku data, i.e. the excess 
variability in the red wing, could be reduced. This would require 
an ionization state of the hot spot under the hard X-ray source 
(where the variable component of the line originates when the 
disc is neutral) corresponding either to resonant trapping of Fe 
photons or to comple te ionization of the Fe atoms (see, e.g., 
l2vcki&Czemvll 19941 for details). 

A remarkable consequence of ionization, crucial for model- 
ing MCG-6-30-15, can be expected also in the soft X-ray band. 
Reflection from ionized material is more efficient at low en- 



ergies than reflection from neutral matter, which may be rele- 
vant in modeling the low energy part of the RMS spectra, typi- 
cally not well described below 3 keV by neutral reflection mod- 
els. Alternatively, changes of a warm absorber could explain 
variability in this band, see next section; distinction between 
these two models with ionized reflection or warm absorption 
is somewhat ambiguous (in soft X-rays), as the photoionized 
disc reflection model reproduces many features in the spectrum 
that could otherwise be interpreted as warm absorption edges 
jCrummv et alj|2006h . 

Our speculative picture for model S (with an ionized in- 
nermost and neutral outer disc) is similar to the model involv- 
ing two distinct reflection regions on the disc, a highly ionized 
inner region at a few Rg and a neutral outer region at larger 
distances, applied to XMM observations of MCG-6-30-15 by 
Ballantvne et al. (2003). Interestingly, such a model can also de- 
scri be the X-ray spectra of oth er high-accretion rate AGNs (see 
e.g. iSobolewska & Donell2007l) . 

The discrepancies predicted by models with vertical mo- 
tion result mostly from the change of the blue peak, which 
should come from a neutral disc (see above). Then, ionization 
effects appear to be not relevant in improving their applicabil- 
ity. Furthermore, in these models irradiation of the disc is much 
more uniform both azimuthally and radially, e.g. for /z, ^ 10 the 
flux irradiating the disc around 6Rg is only larger by a factor of 
2-3 than that around 1 - 2QR„, and, therefore, a range of ioniza- 
tion states along the disc surface, which could give rise to more 
complex variability scenarios, is unlikely. Finally, ionization in 
the inner part, which could be considered to solve t he prob lem of 
the RMS decline around 5.5 keV predicted by the MF04 model 
seems to be ruled out (again due to the very uniform illumina- 
tion) as ionization in that part would lead to strong reduction of 
the whole red wing (which in this model is too weak even for a 
neutral disc, see above), contrary to model S^^, where the az- 
imuthal distribution of the irradiation of the inner disc is very 
non-uniform and a strong red wing can be explained even if its 
variable part is depleted by the strong ionization of a hot spot. 

5.6. Alternative (non GR) models of MCG-6-30-15 

The development of GR variability models was motivated by 
detections of strongly relativistic distortions of the Fe line. A 
further motivation is the fact that the model including a strong 
distant reflection, which could explain some of trends observed 
in MCG-6-30-15, is not applicable to this object (Sect.|42]i. 

The remaining alternative involves a variable absorption 
which can produce an energy dependent RMS variability, with 
larger variations occurring at these energies where the flux is 
affected by stronger absorption. Th e basic picture was devel- 
oped by llnoue & Matsumotol ( l2003h . who postulate changes of 
the column density of an absorber in MCG-6-30-15 on a time 
scale of ~ 10^-10^ sec, under the assumption that the absorption 
cross section does not change. In this scenario, reduced variabil- 
ity between 5 and 7 keV is explained by higher transparency 
of the absorber in this energy range. Similarly, an excess in the 
spectrum in this band - commonly interpreted as a disc line - is 
attributed to the small opacity of the absorber. Note, however, 
that the stronger decrease of the RMS variability around 6 keV 
for larger time bin widths, which motivated the model proposed 
by Inoue & Matsumoto, was seen only in one ASCA observation. 

Several contradictions to such a model were pointed out 
based on the XMM observations (see Vaughan & Fabian 200^ 
for a detailed discussion, disfavouring interpretations of the 
spectral shape excluding the presence of the disc line and rul- 
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ing out explanations of the variability in terms of complex ab- 
sorption). Most importantly, an analysis of the absorption lines 
for the 2001 XMM obser vation did not rev eal any substantial 
changes of the absorber (iTurner et alj |2004|) which appears to 
be constant while the spectrum exhibits changes typical for this 
source. 

Recently, a revised model involving multiple absorption 
zones has been proposed by Miller et al. (2008) and it seems that 
the contribution to the variability coming from changes of the ab- 
sorber remains an uncertain issue. Obviously, photoelectric ab- 
sorption is unimportant above 10 keV and other reduction mech- 
anisms should be taken into account if the reduction of variabil- 
ity ext ends to several tens of ke V, as concluded by Miniutti et^ 
( l2007b . In terms of non-GR models, such a reduction could be 
explained if changes of the spectral index were present with a 
pivot at £■ > 10 keV. Fits to the time resolved spectra with 
the two-component model do not reveal s uch trends (except fo r 
a few cases with the lowest flux, see e.g. iLarsson et al. I l2007h : 
however, they assume the presence of the strong reflection com- 
ponent and, therefore, this conclusion may be model-dependent. 
A model involving short-term variations of the ionization degree 
of warm absorbers, combined with changes of the spectral in- 
dex, for the Suzaku observations is currently under investigation 
(Miyakawa & Ebisawa, in progress). 

5.7. Other objects 

AGNs show a range of spectral variability properties. Some 
show trends similar to MCG-6-30-15 but these are certainly 
not ubiquitous. Accretion rates and black hole masses, similar 
to those in MCG-6-30-15, are typical for Narrow Line Seyfert 
1 (NLSyl) galaxies. A similar geometry of the central region 
may be common in these objects. Indeed, NLSyls often show 
reflection-dominated spectra. Moreover, the blurred reflection 
fits require that emission from within a few Rg dominates, which 
makes them similar to MCG-6-3 0-15; however, in some of 
them the RMS spectra are flat (e.g. IZoghbi et al.ll2008h . On the 
grounds of our results (see Fig.|4^), we point out that such a dif- 
ference may result from slightly different speeds of black hole 
rotation, with the maximal (a > 0.98) speed tentatively expected 
in MCG-6-30-15, and submaximal (a ^ 0.95) which may char- 
acterize other AGNs. An object with the same geometry of the 
innermost region as favoured here for MCG-6-30-15, but har- 
bouring a slightly submaximally rotating black hole, would pro- 
duce a reflection-dominated spectrum, subject to extreme rela- 
tivistic blurring (as r^^ ^ 2), but its RMS spectra would be flat. 

Recent studies of the X-ray spectra in AGNs involving the 
blurred reflection model often indicated rapid rotation of super- 
massive black holes; moreover, the radial emissivity index of the 
fits indicates that emission from the innermost (within a f ew Rg) 
part of the disc dominates (see e.g. ' Crummy et al.ll2006l) . Note 
t hat the relativ istic blurring in such models is typically based on 
a iLaoj (1 19911) line profile, which assumes a maximal value of 
a = 0.998, so the value of a is not measured and the fits re- 
quiring an extreme rotation of the black hole do not exclude a 
slightly submaximal value of a (which, following the above ar- 
guments, may determine the observational diff'erences between 
MCG-6-30-15 and other objects from this class). 

NGC 4051 is an example of NLSyl with a spectral variabil- 
ity simil ar to that in MC G-6-30-15, at least in the higher flux 
states ( Ponti et al.ll2006l) . The two-component model with the 
constant reflection compone nt involving the contribution from 
an ionized inner disc used bv lPonti et al l (L2006) to explain vari- 
ability at higher fluxes is analogous to that describing MCG-6- 



30-15 . Interestingly, in the low flux state studied by iPonti et aTl 
(2006') the variability pattern is different and - in the context of 
GR models - it may imply that the low flux states are related to 
prolonged periods with the hard X-ray source located in the po- 
sition closest to the black hole an d exhibiting mode rate intrinsic 
changes. However, as discussed in lTerashima et al.l (12009). these 
low flux states may be strongly dominated by spectral compo- 
nents formed in distant material, which makes investigation of 
the innermost region less straightforward. 

The RMS spectra decreasing at higher energies, with de- 
pressions around 6 keV , are observed in many AGNs (see e.g. 
iMarkowitz et al.ll2003l) and the two-component model, involv- 
ing the reflection component constant both in normalization and 
in spectral shape, has been popular recently in the modeling of 
such spectra. The first of the observed properties, i.e. the de- 
crease with increasing energy, is predicted by most GR models, 
however, it requires a rather small distance, $ IQRg, of the X- 
ray source, which is not a likely approximation in some cases, 
e.g. for those objects where an extended (~ a few tens of R^), 
hot flow is supposed to form. Moreover, all GR models predict 
some changes in shape and, therefore, results achieved with the 
two-component model cannot be directly interpreted in terms of 
GR effects. In turn, the second of the observed properties, i.e. 
depression around 6 keV, is not easily reproduced by GR mod- 
els and its explanation requires even more extreme parameters, 
namely rapid rotation (with a > 0.95, cf. Fig. 4a) and a distance 
< 3Rg. 

A strictly constant component can be produced by reflection 
from distant matter; however, a strong contributi on of such radia- 
tion is not supported by recent models ( Murphy & Yaqoob 20()9|; 
note that our Fig. [8^ illustrates the strength of effects correspond- 
ing to the amount of reflection from a distant torus estimated in 
that paper, i.e. 6 times weaker than from the pexrav model with 
R - \). Furthermore, such a strong contribution would yield pro- 
nounced, sharp drops at 7 keV both in the RMS spectrum (see 
Fig. ISj)) and in the average energy s pectrum. The latte r is only 
occasionally observed in AGNs (e.g. lBoller et al.ll2002l ; interest- 
ingly, the RMS for this observation is flat). 
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